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Abstract 

We consider a model where a scalar field develops a metastable vacuum state and 
weakly interacts with another scalar field. In this situation we find the probability 
of decay of the false vacuum stimulated by the presence and collisions of particles of 
the second field. The discussed calculation is an illustration of the recently suggested 
thermal approach to treatment of induced semiclassical processes. 



1 Introduction 



The class of problems involving spontaneous and induced semiclassical processes with quan- 
tum fields invariably presents an interesting case study in nonperturbative calculations. The 
most discussed examples of such processes are the Schwinger e'^e" pair creation in an elec- 
tric field [1, 2, 3], the decay of false vacuum [4, 5, 6] and the breakup of metastable strings 
and walls [7]. A semiclassical treatment of such processes can be formulated in terms of 
a configuration in the Euclidean space time, which is a solution to classical field equations 
and is called the bounce [5, 8, 9, 10]. The considered processes can either proceed sponta- 
neously or be induced by the presence of particles and by their collisions, such as in the false 
vacuum decay [11, 13, 14], in the decay of metastable strings and walls [7, 15], and in the 
photon-stimulated Schwinger process [16, 17]. In terms of a Euclidean-space treatment the 
case of induced processes requires one to consider the interaction of the particles with the 
semiclassical configuration. It has been recently suggested that such a consideration can be 
done using an indirect thermal method [7, 18, 19], which appears to be simpler than a direct 
calculation of a relevant correlator [19] in the semiclassical background, especially if more 
than one quantum particle is involved. 

The thermal method allows one to find the rate of the induced process by appropriately 
interpreting the result for the spontaneous process at a small but nonzero temperature, 
where the thermal effects in the rate are contributed by all the n-particle states rather than 
by the vacuum state alone. The contribution of the states with a fixed number of particles 
can be extracted from the thermal rate and thus the decay rate induced by n particles can 
be calculated. An alternative, and more direct, calculation employing correlators of field 
operators in the bounce background is readily tractable in the case of one particle, but 
becomes more cumbersome for processes with a larger number of particles. It has been 
demonstrated [19] that both methods yield the same result for the one photon induced 
Schwinger pair creation. 

In this paper we consider an application of both approaches to the case of interacting 
scalar particles in two dimensions. We consider two interacting scalar fields in the following 
setup. The potential for one of the fields </? is chosen as having the shape shown in Fig.l. 
It is well known that the system with such potential admits quantum tunneling from the 
metastable to the true vacuum, therefore one can discuss the spontaneous decay of the false 
vacuum [4, 5, 8, 6]. Choosing the parameter of the potential asymmetry to be small and the 
mass of the corresponding field </? to be much larger than that of the second field X: M ^ m, 
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greatly simplifies the treatment. Indeed, for a small difference of the vacuum energy density 
e one can employ the so called thin wall approximation [4, 5]. In this case the tunneling 
configuration is given by a disk of the radius 

with the false vacuum outside of the disk and the true vacuum inside. The extent of the 
transition region between the two vacua is assumed to be much shorter than the radius R 
(the thin wall approximation), which is ensured by the condition Mi? ^ 1. The quantity 
ji in Eq.(l) is the mass associated with the transition region, jji ^ M. The smallness of 
the mass m of the second field x compared to M guarantees that in the leading order there 
is no contribution to the interaction of the bubble (bounce) with other objects due to the 
self-action of the field </? but only due to the exchange of the fight particle 

The problem that is addressed here is that of the false vacuum decay stimulated either 
by the presence of one particle of the field x or by a collision between two such particles. 
Previously a similar problem was considered [11, 12, 14] for the particles of the field </?, 
which is always strongly coupled to the bounce in the sense that there are always zero and 
soft modes for this field localized on the bounce. Our present problem is different in that 
we consider the catalysis of the process by particles of a second field x that is only weakly 
coupled to the master field </? so that no anomalously soft modes for the field x arise on the 
bounce. As a result we find closed formulas for the probability of an induced vacuum decay, 
which describe the onset of the expected semiclassical behavior similar to that the catalysis 
by the particles of the master field [11, 12, 14], where the energy of the initial particle(s) is 
transfered to the tunneling degree of freedom. 

Since the main purpose of the present paper is an illustration of an application of two 
approaches to treatment of an induced semiclassical process, we limit ourselves, for simplicity, 
to a model in (1+1) space-time dimensions, while a generalization to higher dimensional 
models is quite straightforward. 

The paper is organized as follows. In Sec. 2 we briefiy describe the model used and 
reformulate the problem in the thin wall limit. In Sec. 3 we find the rate of the metastable 
vacuum decay at small but finite temperature. Using the correspondence between the n- 
particle contribution and the rate at finite temperature we find the width of the x particle 
and the probability of the collision-induced false vacuum decay in Sec.4. In Sec.5 we discuss 
a calculation of the relevant field correlator in the bounce background, and thereby calculate 
anew the decay rate of a single x particle. In Sec. 6 we discuss and summarize the results. 
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2 Vacuum tunneling in Euclidean space time 



We consider a model describing two real scalar fields in 1 + 1 dimensions with the Lagrangian 

^ = ^ id.vf + \ {d,xf - \rn\' - V{^) - M„t(^, x), (2) 

where the potential V^(0) has asymmetric form shown in Fig.l. For further uses this potential 
can be chosen as 

V{^) = \\{v'-v')" -^^{V + V). (3) 

We also choose the interaction potential between the two scalar fields in such a way that it 
is not equal to zero only where the field differs from its vacuum values, for example 

V^nt = a{y - t;^) X = OLp{ip)x , (4) 

with a being a constant. As argued previously, the parameters of the Lagrangian are assumed 
to be such that the mass of the field 93 is much greater then that of the field x, namely 

m<t:M = vV2X. (5) 

In the Euclidean time the solution for the problem of tunneling (metastable vacuum decay) 



Figure 1: The asymmetric potential for the master field (p. 

can be formulated in terms of a semiclassical solution, a bounce. This is a solution to the 
Euclidean equations of motion interpolating between the two vacua separated by a domain 
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wall. The decay rate is found by calculating the partition function around the bounce 
configuration 

I ^ Tf^"^ / ^'^^Xe-'^f^'^l, (6) 
where T, L are the sizes of the system along the Euchdean time and spatial direction respec- 
tively. In the thin wall limit, when the size of the bounce is much greater than the thickness 
of the wall (the inverse mass of the field </?) R — ji/e^ ^/M, one can consider only the soft 
modes of the wall. Then its dynamics is described in purely geometrical terms, namely the 
shape of the wall, which can be represented in polar coordinates by the function r{9) (closed 
trajectory corresponding to the distorted bounce configuration). Therefore the Euclidean 
action of the system can be written in the following form 

■1 ,^ .9 1 



S = fii-eA + J d^x - {d^xf + ^rn\^ + p' 



(7) 



where i — J dOy/r"^ + r^, A — \J dOr"^ are the length of the bounce boundary and the 
area it encircles, and the dot stands for the derivative with respect to 9. The density p^'^^ 
is proportional to the coefficient a and it plays the role of the source for the x field. It is 
obvious that in the discussed approximation the area, where the value of field (p substantially 
differs from the vacuum, is very thin and the support of the function p^'^^ coincides with the 
shape of the boundary of the bounce. Hence, it has the form of a surface 5-function 



p^'>{t,x) = gS{vW+x^ -r{e)), g (x a. (8) 

As a result the partition function can be expressed as an integral over closed trajectories 
describing the shape of the boundary 

I = / '^rVxe-'^''^\ (9) 

with the action S derived from Eq.(7). At a small constant g the latter expression can be 
expanded in a power series in g. Therefore at a small coupling between the scalar fields one 
can find the decay rate of the false vacuum in terms of an expansion in powers of g. The 
decay induced by n particles of the field x arises in the order gr^", while the spontaneous 
vacuum transition proceeds at — > 0. The rate of the spontaneous decay is found from 
calculating the action on the bounce configuration (1) and the determinant corresponding to 
the integration over small fluctuations around the tunneling trajectory, and is given by [6] 
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In the following sections we calculate the rate of the induced decay by using two approaches, 
discussed in the Introduction. 

3 False vacuum decay in thermal bath 

As is already mentioned, one can use a thermal approach to extract [7, 15, 18] from the 
expression for the total process rate at a small but finite temperature the contribution of the 
processes induced by a specific number of particles present in the initial state. Therefore we 
start with calculating the total rate of the false vacuum decay in a thermal bath. It is quite 
clear that at a sufficiently low temperature only the contribution of states containing particles 
of the field x ^i-re relevant, while the contribution from the states with the if excitations are 
suppressed by the Gibbs factor e~^/^. 

At a low temperature one can consider the same effective Euclidean action (7), except 
that now the system lives on the cylinder with the period equal to the inverse temperature 
/3 = 1/T. In other words, one may consider the system living on the (t, x) plane with 
periodic boundary conditions for the Euclidean time coordinate t ^ t + (5 (see Fig. 2). It is 




Figure 2: Periodic copies of the bounce and their interactions through the field x- The self 
interaction (a) within one copy does not depend on the period /3, while the interaction (b) 
between different copies does depend on (3 and describes the thermal effects. 

worth noting that the use of the action from Eq.(7) is justified as long as the temperature 
satisfies the condition T ^ M . However, we consider still much lower temperatures, such 
that a stronger condition T <^ e/fi is satisfied, which drastically simplifies the calculations. 
At such a small temperature (equivalently, at a large period f3) one still may use the same 
bounce configuration to calculate the partition function (9). The zero temperature result 
corresponds to the leading order contribution at /3 — oo and is given by Eq.(lO). The next 
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to leading order corrections arise in fact from the exchange of the x particle, since at T < 2 i? 
there is no bounce deformation [20]. The effect of the exchange of x can be classified into 
two types. One is the self action of the bounce independent of the temperature, it is labeled 
by (a) in the Fig.2. The second, labeled as (b), is the interaction of the two copies of the 
bounce separated along the Euclidean time by /3, 2/3, and so on. In fact, it is still a part of 
the self interaction of one and the same bounce, living on a cylinder, and this part clearly 
depends on the period of the cylinder and disappears when /3 goes to infinity. We restrict 
our interest to the second type of contribution, since this is the only one that produces a 
nontrivial temperature dependence. 

The contribution to the action due to the interaction of two copies of the bounce separated 
by P can be found as usual by using the expression for the propagator of the scalar particle 
with mass m. In 1 + 1 dimensions the propagator D satisfies the equation 

-AL>(i, x) + m^D{t, x) = S^'^\t, x), (11) 

and is given by 

D{t,x) = ^Ko{mr), (12) 

where KQ{mr) is the modified Bessel function of the second kind, and r = + x'^- The 
expression for the part of the action due to the interaction between two copies of the bounce 
then takes the form 

A5'(i) = -i y P^{t\, Xi)D{t2 - h, X2 - Xi)p^{t2 - X2)dtidxi dt2dx2, (13) 
where is the density (8) for the bounce located at the origin 

p^{t,x) = g6{VWT^ - R). (14) 

The shift in the argument of the second function p^{t2 — P,X2) clearly corresponds to the 
second copy of the bounce placed at (/5,0). The extra factor of one half in the expression 
(13) is due to the apparent double counting: the expression without this factor corresponds 
to the interaction between the two copies of the circle, while we need to consider the change 
in the action per one bounce (per one period). 

In order to find the integral in Eq.(13) one can readily notice that the function 

X{t2, X2) = / Do{t2 - h, X2 - xi)p^{ti,xi)dtidxi (15) 
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X{r) 



satisfies the equation 

(-A + m^) x{t, x) = gSiVf^ + x^-R). (16) 

Let (n, 6i) be the polar coordinates corresponding to the Cartesian ones {U, Xi). The solution 
to the equation (16) can then be found as 

gRIo{mR) Ko{mr), r > R, ^^^^ 
gRlQimr) KQ^mR), r < R, 

and the integration over the first circle then reduces to the following 

J Ko (m^r^ + R'^ - 2r2R cos dOi = 27r Io(mR) Ko(mr2) , (18) 

where Io{mR) is the modified Bessel function of the first kind. Performing the further 
integration over r2 and 62 in Eq.(13) and summing over the contributions from the interaction 
with all copies one finally finds 



Ft e 
— = — exp 
L 27r ^ 



TTfJ, 



AS 



(19) 



with the total temperature dependent part of the action given by the following expression 

00 

AS = -2Trg^R^I^{mR) ^ Ko{ml3n), (20) 



n=l 



Expanding the expression (19) in AS at large /3 one can find the low temperature behavior 
of the thermal corrections to the decay rate: 



L 



27r 



exp 



l + 2ng'R^mmR) 



TT 



2m 13 



e-'^P + . . . 



(21) 



As expected, these corrections are exponentially suppressed in the parameter m/T. 



4 Induced false vacuum decay 

The enhancement of the tunneling at a finite temperature arises through the stimulation 
of the process by the particles present in the bath. The dependence of the rate of the 
process induced by the x particles on their energies then translates into the dependence on 
the temperature T after averaging over the thermal distribution of the particles with the 
standard density function 
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where the energy is given by — y k"^ + m?- The number of particles involved in each of 
the microscopic processes can be readily identified by the power of the factor g^. Since the 
thermal correction (20) in the action is proportional to g^, the n-particle contribution to the 
thermal rate is given by the n-th power of AS* in the expansion of the factor exp (—AS") in 
the expression (19). In particular, the one and two-particle contributions to the decay rate 
in a thermal state are given by 



2X 



e 


exp 


Tlfl 


2^ 


e 


e 


exp 


nfj, 


2^ 


e 



(23) 



On the other hand, the one-particle contribution can be found in terms of the width of 
the particle x associated with the destruction of the false vacuum by the particle's presence: 



+°° dk 



27r ^Jk'^ -\- rm? 



^ n=l -^"^ 



dx 



e — 



TT 



Y,Ko{mi3n). 

n=l 

(24) 

The origin of the Lorentz factor m/\Jk? -\- m? is rather straightforward: the width 7^ is 
measured in its rest frame, while we consider the process in the laboratory frame. Comparing 
the two results (23) and (24) we find the width 

,2" 



7x 



n e 



m 



g^WI^{mR) exp 



TT/X 



(25) 



Similarly, one can express the two-particle contribution from (23) through the rate 
w{ki,k2) of the metastable vacuum destruction by collisions of two particles. Namely, the 
relation is as follows 

dkidk2 



(26) 



The two-particle rate function w is related to the dimensionless and Lorentz invariant 1-1-1 
dimensional analog cr of a cross section in 3 -|- 1-dimensions: 



w{ki, k2) = cr{ki, k2) x Flux = cr , 



(27) 



where Vrei is the relative velocity of the two particles, which is expressed in terms of the 
energies of the particles and of the standard kinematical invariant / = \/ {ki ■ /c2)^ — fn^ as 

/ 



Vrel = 



(28) 
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Due to the Lorentz invariance, the function a{ki, /C2) can depend only on the invariant /, 
so that one can consider it in a form of a power series in /: 



(29) 



In order to write the integral over the momenta of the two particles with proper inte- 
gration limits one should take into account that the particles are identical and also that for 
the collision to take place they should move toward each other. Let particle '1' be the one 
on the left and the particle '2' on the right. Also let the spatial momentum k be considered 
as positive when the particle moves from left to right. Then the coUision takes place if and 
only if ki > k2, where both the spatial momenta can be of either sign, and the integral in 
Eq.(26) takes the form 



Ik ^ r 

L i-, 



dki 
2nE, 



n 



/ki 
-00 



dk2 



27: El 



k2 



n{k2)a{I)I . 



(30) 



This integral can reproduce the two-particle thermal term in Eq.(23) only if the product <y{I)I 
is constant, i.e. only if the expansion in Eq.(29) reduces to a single term (t(/) = c_i//. The 
integrand in Eq.(30) then factorizes into terms each depending on the absolute value of the 
momentum, so that the limits of integration can be rearranged as 



2x 



C-1 



/+00 
-00 



dki 
27: E, 



n{ki) 



dko 



dki 



27: Ei 



ki 



n{ki 



00 27r £'fe2 
dk2 



n{k2) 



27: E, 



k2 



n{k2 



2c_ 



E^o(m/3n) 



.n=l 



(31) 



Upon a comparison with the expressions (20) and (23) one thus arrives at the final result for 
the dimensionless invariant 'cross section' for the destruction of the false vacuum in coUision 
of two x-particles: 



a{I) 



7re 
21 



27:g^R^I^{mK) 



exp 



TT/i 



(32) 



5 Direct semiclassical calculation for one particle 

In this section wc consider another approach to calculating the rate of the false vacuum decay 
induced by one particle, which is similar to the one used in Rcf. [19] for the photon-induced 
Schwinger process. Interpreting the particle-induced decay rate as the width of the particle 
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X, one can find it from tlie propagator of the x fi^ld in the background of the bounce as the 
imaginary part of the density- density correlator^: 

n(t,a;) = (p(t,x)p(0,0)). (33) 

Using the formulas (7) and (9) we can rewrite the leading semiclassical term of the correlator 
in the following form 

Imn(t,a;) = Im ^ r'^p(t, x)p(0, 0) e"^['^l j ^^^P^'"H^, a;)p(^')(0, 0)e-^M 

= JJ^ J P^it-tB,X - XB)p^i-tB,-XB)d'^XB, (34) 

with p^{t,x) from Eq.(14). 

The integration in Eq(34) is performed in the following manner. We choose the coordinate 
system such that the coordinate x in the integral is equal to zero. One can always do so by 
an appropriate rotation of the coordinates. Due to the 0(2) symmetry, the full final result 
can be restored in an arbitrary coordinate system by simply promoting t to y/W+lc^ in the 
final expression for the correlator. After this choice of orientation we make the change of 




Figure 3: The two positions of the bounce for which the density correlator is not equal to 
zero. 

the variables of integration 

zi = \/tl + xl, 

Z2 = ^{tB-ty + xl, (35) 

^This is in a one to one correspondence with QED, where the attenuation rate of the photon is given by 
the imaginary part of current-current correlator. 
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and rewrite the integral in the following form 

/ J dtBdxBS (y^tl + xl - 5 (^^{ts-tf + xl ~ R^ 

= 2g^ [ dz,dz25{z, - R)S{Z2 - R)J f ^) dz,dz2 = f^'f . (36) 
J \ Z\Z2 / ty AR"^ — t 

where J (^^ff) is the Jacobian of the transformation (35). The factor of two in the second 
expression in Eq.(36) comes from two possible configurations shown in Fig. 3. After promoting 
t — > + a;^ and making the Fourier transformation one finds 

where the two-momentum q is set on the mass shell, = —w?. Using the formula (34) one 
thus finds: 

(38) 

which coincides with Eq.(25). 



7x = = — 9 R {mR) exp 

mm 



6 Discussion and summary 



The formulas for the false vacuum decay rate induced by one (Eq.(25) or Eq.(38)) and 
by a collision of two (Eq.(32)) light particles of the weakly interacting field x merit some 
remarks. Firstly, one can readily notice that the width of the particle associated with the 
destruction of the vacuum is singular in the limit m ^ 0, which is a consequence of the 
infrared behavior in this limit in 1+1 dimensions. The singular growth at low m however 
does not result in any unphysical artifacts, since the rate of decay for a state with any fixed 
finite energy E is finite in the massless limit due to the time dilation factor m/E. Secondly, 
at large mR the exponential asymptotic behavior of the Bessel function correctly reproduces 
the known [11, 12] semiclassical exponential factor in both the one-particle and two-particle 
processes: 



7x ^ exp 



+ 2mR 



w ~ exp 



TTyU 



AmR 



(39) 



It is worth noting that the enhancement factor in the collision is determined by the mass of 
the particles rather than by their cm. energy ^/s. This behavior, which is in agreement with 
the geometrical calculation of the exponent in Refs. [11, 12] is a consequence of the assumed 
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linear in x interaction in Eq.(4). Had we assumed instead a quadratic in x interaction, the 
exponential factor would be determined by -^/s R at -/s ^ m, as is the case for the similar 
factor in the rate of destruction of metastable strings and walls by Goldstone bosons [7]. 
Moreover, the absence of the dependence on the energy of the exponent in a is true only 
as long as the energy is small in comparison with the mass scale M of the master field (/?, 
^/s <C M . The reason for this behavior is that iteration of the interaction in Eq. (4) through 
the short-distance propagator of (/? gives rise to terms with higher powers of x, which however 
are suppressed by inverse powers of M . 

In summary. We have considered an illustrative 1+1 dimensional model of interacting 
scalar fields and demonstrated an application of the semiclassical thermal method to calcu- 
lation of the rate of destruction of the false vacuum by one- and two- particle states of a light 
scalar field. For the one-particle state a direct calculation of the rate is also done in terms 
of the imaginary part of the particle propagator in the tunneling background. In the latter 
case the direct method is possibly simpler, while in the case of two particles the thermal 
approach appears to be more straightforward. Furthermore, the formula for the thermal 
decay rate produces a generating function for the probability of the process induced by an 
arbitrary number of particles. 

The calculations discussed in the present paper can be readily generalized to other models 
of interaction and to models in higher dimensions. 
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